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quark actions. 
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1. Introduction 

The decay processes of the neutral kaon system contain rich physics in the SM and beyond 
it. The indirect CP violation parameter e provides us a chance to examine the SM and make 
constraints on the physics beyond it. In order to determine e we need the kaon matrix element of the 
effective Hamiltonian for AS = 2 transition, whose general structure is written as (K°\H^ =2 \K°), 
where H^= 2 = G F /\/2^ ; - V^C^/l)^- with G F the Fermi constant, the CKM factors, C, 

the Wilson coefficients and G\ the relevant local operators for the process. Although it is known 
that only the left-left operator is relevant in the SM, we need the operators with more general 
chiral structures for the physics beyond the SM. For example the supersymmetric model requires^] 

0\ = s a Y^-Y5)d a s b Y^l- Y5 )d h , (1.1) 

G 2 = s^(l-y 5 )d a i } {l-y 5 )d\ (1.2) 

^3 =^(1-75)^(1 -75K, (1-3) 

^4 = ^(1-75)^(1 + 75)^, (1.4) 

^5 = 5°(1 - 75 M V ( 1 + 75) d a (1.5) 

and ^1,2,3 obtained from the ^1,2,3 by exchanging the left- and right-handed quarks. 6\ is the same 
operator in the SM, while ^2,3,4,5 an d G\ .2,3 are the non-SM operators. 

In this report we present the results of (K°\ €?j\K } (i= 1 , . . . , 5) obtained in quenched domain 
wall QCD. Although these operators are introduced by the supersymmetric model, we should note 
that (K°\ &i\K } with i = 1, ... ,5 are a complete set of the kaon matrix elements of the four-quark 
operators allowed by the symmetries. This implies that they are sufficient for analysis of any model 
beyond the SM. We are particularly interested in the relative magnitude of non-SM to SM matrix 
elements: 

R(u)- <*V'fr)l*°> i- 2 5 H6) 
* ,(M)= <KO|*0z)|*o>' i - 2 "-" 5 ' (L6) 

We also present the kaon B-parameter for the SM: 

B lW = B KW = 8/g0 ,_ ,, A , /A |_ "777^7' ( L7 ) 

and those for the non-SM: 

^^ J -Q(^ k-75rf(M)|0)(0|^(M)l^ )' ' U ' 

where Q = {|,— i,— 2, — §} are the convention factors introduced by the vacuum saturation ap- 
proximation. Our results are compared with the previous works employing the overlap [||] and the 
clover quark actions^]. 

2. Simulation details 

2.1 Simulation parameters 

In Table |l| we summarize the simulation parameters, part of which are the same as in the 
previous CP-PACS calculation of £'/ £ lFll- The domain wall fermion with the wall height M = 1 .8 is 
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Table 1: Simulation parameters. P and R denote the expectation values for the plaquette and the 1x2 
rectangular. 





plaquette 


Iwasaki 


p 


6.0 


2.6 


N s xN t xN 5 


16 3 x 32 x 16 


16 3 x 32 x 16 24 3 x 32 x 16 


a _1 [GeV] 


2.12@| 


2.00[|] 


P 


0.59374[@] 


O.67O632(1O)[0] 


R 




0.45283(2)|§ 


No. config. 


180 


400 200 






(400 for m f a = 0.02) 


m 5<? [MeV] 


3.14(36) 


0.292(44) 



employed with the plaquette and the Iwasaki gauge actions in quenched approximation. We choose 
N s x N t x N5 = 16 3 x 32 x 16 lattice at /3 = 6.0 for the plaquette gauge action, while N„ = 16 3 
and 24 3 at jS = 2.6 for the Iwasaki gauge action in order to investigate the volume dependence. 
As for the lattice spacing we refer to Refs. [j5L |6|] where the lattice spacing determined from the 
Sommer scale ro = 0.5fm is parameterized in terms of j3. We take degenerate masses rrifa = 
0.02,0.03,0.04,0.05,0.06 for the up, down and strange quarks in both gauge actions. For each 
value of mfa gauge configurations are independently generated with the algorithm of the five-hit 
pseudo-heat-bath combined with four overrelaxation sweeps, which we call an iteration. In Table 
we list the numbers of configurations with the interval of 200 iterations for measurements. The 
anomalous quark masses m$ q for both gauge actions are calculated in Ref.[||]. Their values in the 
physical unit given in Table are obtained by using the lattice spacing discussed above. 

2.2 Calculational method 

The ratios R{ are extracted from the plateau of 

choosing < f < iV ( - 1. Similarly Bk and fi, (i = 2, • • • ,5) are obtained from the three-point 
functions divided by two-point functions, 

^(K°(z,N t - m(y,t)K°(x,0)) 
lLly(K°(z,N t - l)A {yj))Z f f(A (y',t)K (x,0)) 
^ sr AK°(z,N t -l)^(y,t)K (x,0)) 

QZ?$(K (z,N t - i)P(y,t))'LfA p (y , ,t)K (x,o)} 

In order to calculate the hadron correlation functions we solve quark propagators with the Coulomb 
gauge fixing employing wall sources placed at the edges of lattice where the Dirichlet boundary 
condition is imposed in the time direction. 

We estimate errors by the single elimination jackknife procedure for all measured quantities. 



for B K , 



(2.2) 



for B:. 



(2.3) 
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2.3 Renormalization factors 

In order to convert the matrix elements obtained on the lattice to those defined in the contin- 
uum, we make a one-loop perturbative matching of the lattice operators and the continuum ones 
at a scale \i = \ja, where the latter is defined in the MS scheme with the naive dimensional reg- 
ularization(NDR). Since the domain wall quarks retains the good chiral symmetry on the lattice, 
the four-quark operators are renormalized in the same way as in the continuum: multiplicative 
renormalization for G\, while mixing between &2 and ^3 and between and 

The one-loop perturbative renormalization factors with mean field improvement are written 

asfBl 



,lat 



(2.4) 



where wo is a function of domain wall height, Z w is the renormalization factor of domain wall height 
and Zij are those of the four-quark operators. Indices i, j label the operator number. Numerical 
values of Zn with mean field improvement for the Iwasaki and the plaquette gauge actions are 
given in Ref.[10]. Employing j3 = 6.0, cT x = 2.12GeV and P = 0.59374 for the plaquette gauge 
action we replace M — > M = M + A(u - 1) = 1.311 with u = P 1 / 4 and obtain 



/ 0.7287 
0.6845 



^plaquette 



(ji = 1 /a) 



\ 









0.0352 







-0.00156 
0.8682 











0.6325 







-0.0414 



(2.5) 



-0.0689 0.7564 / 



For the Iwasaki gauge action we obtain M = 1.420 and 



/ 0.8062 
0.8124 



^Iwasaki 

"'7 



(M = 1/a) 



V 









0.0156 







-0.00679 
0.9241 











0.7847 







-0.0427 



(2.6) 



-0.0477 0.8425 / 



for p = 2.6, a- 1 = 2.00GeV, P = 0.670632(10) and R = 0.45283(2). The details of the mean 
field improvement for the domain wall fermion are explained in Refs.[^, [T(J. We observe that 
the Iwasaki gauge action shows smaller one-loop corrections than the plaquette gauge action. For 
the axial vector current and the pseudoscalar density in Bk and Bj (i = 2, . . . ,5), we also use the 
one-loop perturbative renormalization factors with mean field improvement. 



3. Results for R; and B t 



Since the kaon matrix element of 0\ is proportional to M\, the ratios Rj in Eq.(|L6|) should 
diverge toward the chiral limit. In order to keep the ratios finite in the chiral limit and tame the 
quark mass dependences of /?,, we introduce another definition: 



1 



(K°\^)\K°) 



(3.1) 



lat 
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Table 2: Results for Ri and B, extrapolated at the physical point with the quadratic fit. The renormalization 
scale is fi = I /a 



N s xN t x N 5 


plaquette 
16 3 x 32 x 16 


Iwasaki 

16 3 x 32 x 16 24 3 x 32 x 16 


Ri 
R 2 
R 3 
R4 
R5 


1 

-18.22(69) (38) 
4.74(17)(10) 
27. 19(82) (24) 
8.11(24)(6) 


1 

-19.65(45)(21) 
5.21(12)(6) 
30.14(55)(11) 
9.01(16)(3) 


1 

-18.97(16)(16) 
5.039(41)(44) 
29.68(20)(13) 
8.794(58) (34) 


Bi 

B 2 
B 3 

B 4 
B 5 


0.561(13) 

0.5603(80)(9) 

0.731(H)(1) 

0.691(13)(8) 

0.616(12)(7) 


0.54361(65) 
0.5564(52)(10) 
0.7368(71)(15) 
0.7075 (72) (31) 
0.6346(65)(25) 


0.5233(55) 
0.5369(17)(3) 
0.7124(24) (4) 
0.6996(22) (26) 
0.6226(19)(24) 



where dim is the pseudoscalar meson mass on the lattice and M^ xp is the experimental kaon mass. 

In Figure [j] we plot the results for Ri as a function of m 2 M . We observe little volume dependence 
for the Iwasaki gauge action. The results of the Iwasaki and the plaquette gauge actions show 10% 
discrepancies, which could be 0(a 2 ) effects due to the difference of the gauge actions. We extrap- 
olate the results for Ri at the physical kaon mass, which is located around mja ~ 0.02, employing 
the quadratic and the logarithmic fitting functions: 

Ri = b + b l m M + b 2 (m M ) 2 , (3.2) 
Ri = b + b\m M + b 2 m M \og(m M ). (3.3) 

The results are depicted in Figure [T], where both fitting functions give consistent results at the 
physical point. Table ^ summarizes the values of Ri extrapolated at the physical point with the 
quadratic fitting function. The first errors are statistical. The second denote differences of the 
central values for the two fitting functions, which are considered to be the systematic errors. We 
check numerically that the uncertainties due to a choice of the fitting functions are smaller than 
the statistical errors. The results for Rj show a striking feature that the non-SM matrix elements 
are much larger than the SM ones for both gauge actions. Especially, the magnitudes of the matrix 
elements {K°\^ 2 (n)\K°) and (K \<ff 4 (li)\K Q ) are one order larger than that of {K°\0i(p)\K°). 

In Figure || we also show the results for the B-parameters denned in Eqs.(L7) and ( |1.8[ ). Both 
gauge actions give consistent results ranging from ntf-a = 0.02 to 0.06. We again find little volume 
dependence for the Iwasaki gauge action. In order to extrapolate Bk at the physical kaon mass, we 
use the following fitting function suggested by ChPT: 

B K = B(l-3cm M \og{m M )+bm M ). (3.4) 

As for Bi (i = 2, . . . , 5) we employ two fitting functions: 

Bi = b + bim M + b 2 (m M ) 2 , (3.5) 
Bi = b + bimM + b 2 m M log(m M ). (3.6) 
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Figure 1: Rj as a function of m 2 M . Solid symbols denote the extrapolated values at M^ xp . The lattice spacing 



is fixed by the Sommer scale 



Figure || shows that the extrapolated values at the physical point with both fitting functions are 
consistent. This is also confirmed by the numerical values listed in Table ||, where the first errors 
are statistical and the second for differences of the central values for the two fitting functions. 



4. Comparison with previous works 

In Figure [j] we compare our results for Rt with those obtained using the overlap fermion and 
the plaquette gauge action^]. The lattice spacing is fixed by the Sommer scale ro = 0.5fm. Both 
results show good consistencies for R2 and Rj, while large deviations are observed for /?4and R$. 

As for the B-parameters we plot in Figure |2| the previous results with the overlap fermion[||] 
and the 0(a)-improved Wilson quark action^]]. We again employ the Sommer scale to determine 
the lattice spacing. Large discrepancies are observed for all the B-parameters except B$, though the 
previous results have rather large statistical errors. 

Although the different choices of the quark and the gauge actions should allow the 0(a 2 ) un- 
certainties, the magnitudes of the inconsistencies between the previous results and ours are more 
than expected. We suspect that a possible source of discrepancies is difference of the renormal- 
ization methods: one-loop perturbation for our results and nonperturbative MOM scheme for the 
previous ones. We are now working on a nonperturbative renormalization for the domain wall 



fermion using the Schrodinger functional method [|1 1 
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Figure 2: Bag parameters B, as a function of m M . Solid symbols denote the extrapolated values at M K . The 
lattice spacing is fixed by the Sommer scale. 
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